Exploration on the relativistic symmetry by similarity renormalization group 
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The similarity renormalization group is used to transform Dirac Hamiltonian into a diagonal form, 
which the upper(lower) diagonal element becomes an operator describing Dirac (anti)particle. The 
eigenvalues of the operator are checked in good agreement with that of the original Hamiltonian. 
Furthermore, the pseudospin symmetry is investigated. It is shown that the pseudospin splittings 
appearing in the non-relativistic limit are reduced by the contributions from these terms relating the 
spin-orbit interactions, added by those relating the dynamical terms, and the quality of pseudospin 
symmetry origins mainly from the competition of the dynamical effects and the spin-orbit interac- 
tions. The spin symmetry of antiparticle spectrum is well reproduced in the present calculations. 
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Many years ago a quasidegeneracy was observed in heavy nuclei between single-nucleon doublets with quantum 
numbers (n,l,j = I + 1/2) and (n — 1,1 + 2,j = I + 3/2) where n, I, and j are the radial, the orbital, and the 
total angular momentum quantum numbers, respectively P, Q. The quasidegenerate states were suggested to be 
pseudospin doublets j = I ± s with the pseudo orbital angular momentum I = I + 1, and the pseudospin angular 
momentum s — 1/2, and have explained a number of phenomena in nuclear structure. Because of these successes, 
there have been comprehensive efforts to understand the origin of this symmetry. Until 1997, it was identified as a 
relativistic symmetry Nevertheless, there is still a large amount of attention on this symmetry. The pseudospin 
symmetry (PSS) of nuclear wave functions was tested in Refs.0, H| with conclusion supporting the claim in Ref.Q. 
The existence of broken PSS was checked in Rcfs.[6, 7], where the quasidegenerate pseudospin doublets were confirmed 
to exist near the Fermi surface for spherical and deformed nuclei. The isospin dependence of PSS was investigated 
in Ref. where it is found that PSS is better for exotic nuclei with a highly diffuse potential. PSS was shown to 
be approximately conserved in medium-energy nucleon scattering from even-even nuclei [9h11|. In combination with 
the analytic continuation method, the resonant states were exposed to hold the PSS in Refs.[12l. Il3|. In Ref. 14 1, the 
conditions which originate the spin and pseudospin symmetries in the Dirac equation were shown to be the same that 
produce equivalent energy spectra of relativistic spin- 1/2 and spin-0 particles in the presence of vector and scalar 
potentials. Furthermore, the symmetries and super-symmetries of the Dirac Hamiltonian were checked for particle 
moving in the spherical or axially-deformed scalar and vector potentials [15j ]. More reviews on the PSS can be found 
in the literature [l6| and the references therein. Recently, a perturbation method was adopted to investigate the spin 
and pseudospin symmetries by dividing the Dirac Hamiltonian into the part of possessing the exact (pseudo)spin 
symmetry and that of breaking the symmetry 17 1 . 

Despite the large number of studies on PSS, it is still not fully understood the origin of PSS and its breaking mech- 
anism since there is no bound states in the PSS limit. Hence, many efforts are devoted to compare the contributions 
of different terms in the Schrddinger-like equation for the lower component of Dirac spinor to the pseudospin energy 
splitting. In Refs. [l8l - [20l | . the PSS in real nuclei was shown in connection with the competition between the pseudo- 
centrifugal barrier and the pseudospin-orbital potential. In Refs. 21, 2^, it was shown that the observed pseudospin 
splitting arises from a cancellation of the several energy components, and the PSS in nuclei has a dynamical character. 
A similar conclusion was reached in Refs. 23j, [24j. However, in these studies, one encounters inevitably the singularity 
in calculating the contribution of every component to the pseudospin splitting, and the coupling between the energy 
e and the operator in solving the Schrddinger-like equation for the lower component of Dirac spinor (to see Eq.(j4]) 
in the following), which affect our understanding on the origin of the PSS. As seen in the following Eq.(|4]), it seems 
that only f^jp- destroys the PSS, but the pseudospin splittings are related to every component 21-24j]. In order to 
cure these defects, in the paper we transform the Dirac operator into a diagonal form by similarity renormalization 
group (SRG), in which the upper(lower) diagonal part becomes an operator describing Dirac (anti)particle with the 
singularity and the coupling disappearing. In the following, we first derive out the operator, and then present its 
application in analyzing the PSS. 

Assuming the spherical symmetry, the radial Dirac equation can be cast in the form of 



H s ip = op, 
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where E(r) = V(r) + S(r) and A(r) = V(r) — S(r) denote the combinations of the scalar potential S(r) and the 
vector potential V(r), and K is defined as K = (I — j)(2j + 1). To understand the PSS, one decouples Eq.© into two 
equations for the upper and lower components: 
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here the effective masses 2M + = e + M — A and 2M_ = e — M — S. The prime denotes derivative with respect to r. 
From Eq.Q, it can be seen that the system possesses exact PSS when £' = 0. Unfortunately, the condition cannot 
be realized in real nuclei, many efforts are devoted to analyze the contributions of various terms to the PSS [2l|-[24j]. 
However, as there exist deficiencies mentioned before, we decouple Eq.([T|) by SRG. 

Without loss of generality, we begin our formalism for a general Dirac Hamiltonian H = a -p + j3(M + S) + V, which 
is fully applicable for H s . Following Wegner's formulation of the SRG 25], the initial Hamiltonian H is transformed 
by the unitary operator U (I) according to 



where I is 



H (I) = U (I) HU^ (1) , H(0)=H 
a flow parameter. Differentiation Eq.(|5]) gives the flow equation as 



dl 



H(l) = [r,(l),H(l)], 



(5) 
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with the generator r/(l) = dU ^ (I). There are several possibilities to choose the T}{1) so that H (I) becomes diagonal 
in the limit I oo. In Wegner's original formulation 25j, rj(l) was chosen as the commutator of the diagonal part of 
H{1) with H(l) itself, i.e. r/(l) = [H^i ag (l), H(l)]. An alternative to Wegner's formulation is r/(l) = [G,H(l)], where G 
is a fixed (/-independent) hermitian operator. It is straightforward to show that H(l) converges to a final Hamiltonian 
which commutes with G. Here, we hope to transform Dirac Hamiltonian into a diagonal form, which must commute 
with the (3 matrix. Thus, it is appropriate to choose rj(l) in the form 



V (l) = [PM,H(l)}. 



(7) 



In order to solve Eq.©, the technique in Ref.[26( is adopted. The Hamiltonian H(l) is presented as a sum of an even 
operator s(l) and odd operator o(l): 



H(l) = e(l) + o(l), 



(8) 



where the even or oddness is defined by the commutation relations of the respective operators, i.e., e(l)/3 = (3e(l) and 
o(l)fi = -/3o(l). To put Eqs.© and © into Eq.© gives 
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The system of Eqs.|9]) and ([TOf can be solved perturbatively in 1/M [26]. It is convenient to introduce a dimensionless 
flow parameter A = IM 2 . Since e(0) = j3 (M + S) + V, and o(0) = a ■ p and the expansion of e(X)/M in a series in 
1 /M contains terms starting with the zeroth order term 
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whereas the expansion of o(A) /M starts with the first order 

1 00 1 



J'=l 



Differentiation Eqs. ([TT|) and (|f 2p yield the following equations, 
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The solutions of the equations (|T3[) and (1141 are obtained as 
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with the initial conditions 



e„(0) = /3,£i(0)= ( 85 + V,£ n (0) = if n ^ 2, 
01 (0) = a-p,o n (0) = if n> 2. 



(17) 



From the equations (|15H16|) with the initial condition (TT7)) . we obtain £o(A) = /?, £1 (A) = /3S + V, and 01 (A) = 
o\ (0) e _4A . Hence, it is easy to verify that o„(A) exponentially goes to zero when A — > 00. So, the diagonalized Dirac 
operator is obtained as 

e(oo) = Me (oo) + £1(00) + -^£ 2 (oc) + -^£3(00) + -^£ 4 (oc) H 

= M£ (0) + £! (0) + 2^/80? (0) + g^j [[Ol (0) , £l (0)] , 01 (0)] 

+ (-40 4 (0) + 01 (0) [[Oi (0) , £! (0)] , £l (0)] + [[ 01 (0) , £1 (0)] , £l (0)] 0l (0) - 2 [ 01 (0) , £1 (0)] [ 01 (0) , £1 (0)]) + 

Here, only a spherical system is considered, £i(0) = 
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is an operator describing Dirac particle with p 2 
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is an operator describing Dirac antiparticle with p 2 = — + ^ • The Hamiltonian for Dirac antiparticle in this 
case is in fact —H2 with eigenvalue — e. This is consistent with the transformations of potentials S — > S, £ — > — A, 



and A — ► — £ under charge conjugation from Hi 22j. Here, the primes have the same meaning as that in Eq.(j4]) and 
the double primes denote second-order derivatives with respect to r. 

The first two terms of Hi{H2) correspond to the operator describing Dirac (anti)particle in the non-relativistic 
limit. The relativistic effect begins to show up from the order of 1/M 2 in the perturbation expansion of s(oo), which 
are presented in H\(H2) from the third to fifth terms. In order to obtain better result, the perturbation expansion 
up to order 1/M 3 is also included in Eqs. (fl"9|) and ([20]) . 

In Eq. fTi?)) , it can be seen the spin symmetry is exact for Dirac particle when A' = 0, which agrees with Eq.©. 



The same result is obtained for Dirac antiparticle in comparing with Ref. [22L 1271] . which can be observed from Eq. ([2H|) 
with £' = 0. Particularly, the singularity disappears in every component of Eqs. (fT9|) and (f20|). the operators Hi and 
H2 are Hermitian. In addition, there is no the coupling between the energy e and the operator Hi{H2). Thus, the 
energy spectra of H\ and H2 can be calculated conveniently. 

The energy spectra of Hi(H 2 ) agree the results of Eq.([T]) very well, and the energy splittings of pseudospin partners 
are in agreement with the exact relativistic case. Especially, the contribution of every component to the pseudospin 
splittings can be calculated, which is helpful to analyze the origin of PSS. In order to convince the conclusion, a 
Woods-Saxon type potential is adopted for £(r) and A(r), i.e., £(r) = £o/(as, rs, r) and A(r) = A /(aA, »"a, r ) 
with 

/(a , r , r) = — r-. (21) 
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The corresponding parameters are determined by fitting the energy spectrum from the RMF calculations for 208 Pb 
(to see Ref. [HI)- The energy spectra of Hi(H2) are calculated by expansion in harmonic oscillator basis. The energy 
spectra of the six pseudospin partners are shown in Fig.l, where the first column in each subfigure corresponds to 
that H\ is approximated to the non-relativistic limit. The second and third columns in each subfigure correspond to 
that Hi is approximated to the order 1/M 2 and 1/M 3 , respectively. The exact relativistic spectra (the eigenvalues of 
Eq-(HI) are displayed in the fourth column. From Fig.l, it can be seen that the deviations between the non-relativistic 
limit (the first column) and the exact relativistic case (the fourth column) are very large, i.e., the relativistic effect is 
apparent in the present system. With the increasing perturbation order, the calculated result is closer to the exact 
relativistic one. When Hi is approximated to the order 1/M 3 , the calculated spectra are considerably agreeable with 
those from the exact relativistic calculations. Especially, the pseudospin energy splitting is in good agreement with 
the exact relativistic result. These show the operator Hi presents a good description for Dirac particle, and can be 
used to analyze the PSS. 

In order to disclose the origin of PSS, we check the contribution of every component in Hi to the pseudospin energy 
splitting. Based on the consideration of attribute and hermitian, we decompose Hi into the eight components: £ (r) + 

6. -5jfr (V - S '£)> 2^ (V " 2S'f), - S ' 2 - 2 ^ 4SS " > "A, which are respectively 

labelled as Oi, 0%, ■ ■■ ,0$. For the fc-state with eigenvector ipk, the contribution of Oi to the level Eu is calculated 
by the formula (k\ Oi \k) = ifjlOiipkd 3 f 7 which is denoted as ei(k). To reduce the length of the article, only the 
data for two pseudospin doublets are listed in Table I. 

The contribution of Oi to Ae is very large, which means that the spectra of Hi in the non-relativistic limit do 
not have the PSS. Namely, PSS is not a non-relativistic symmetry. It has a relativistic origin, which agrees with 
the claim in Ref.Q. The contributions of O3 and Oq to Ae are negative, which implies the pseudospin splittings in 
the non-relativistic limit are reduced by the contributions from these terms relating the spin-orbit interactions, and 
agrees with the relativistic interpretation of PSS. The contributions of O2 and O5 to Ae are positive, which means the 
pseudospin splittings are added by the contributions relating the dynamical terms, and supports a relativistic origin 
for this symmetry. Compared with O2, O3, O5 and Oq, the contributions of O4, Oj and Og to Ae are relatively minor. 
These show the quality of PSS origins mainly from the competition of the spin-orbit interactions and the dynamical 
effects. Although the contributions of O4, O7 and 0% to Ae are minor, their influences on PSS can not be ignored, 
which supports partly with the claim in Refs.[2l|-[23j], i.e., the observed pseudospin splitting arises from a cancellation 



of the several energy components, and the PSS in nuclei has a dynamical character. 

In order to check further the applicability and validity for the present formalism, we have calculated the energy 
spectrum of H2 for Dirac antiparticle, which is shown in Fig. 2. The good spin symmetry is displayed clearly, which 



is in agreement with Refs. [22L l27j. 



In summary, the similarity renormalization group is used to transform the spherical Dirac operator into a diagonal 
form. The upper(lower) diagonal element becomes an operator describing Dirac (anti)particle, which holds the form of 
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TABLE I: The contribution of the operator Oi to the level Ek'- e»(fe) = {k\ Oi \k), where Oi can be seen in text, k represents 
the single particle states 2d 5 /2, 137/2, 2/7/2, and lhg/z, and Ae = ei(a) — The data listed in the last line is a sum from 

the first line to the eighth line. 
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FIG. 1: (Color online) The energy spectrum of Hi for the six pseudospin partners. The first column in each subfigure corresponds 
to that Hi is approximated to the non-relativistic limit. The second and third columns in each subfigure correspond to that 
Hi is approximated to the order 1/M 2 and 1/M 3 , respectively. 



Schrodinger-like operator with the singularity disappearing in every component. The energy spectra of the operator 
are calculated in good agreement with the exact relativistic ones. By comparing the contributions of the various 
components to the energy splittings, PSS is shown to be a relativistic symmetry. The quality of PSS is correlated 
with the contribution of every component of Hi to the pseudospin splitting, especially, the competition of the spin- 
orbit interactions and the dynamical effects, which supports the claim of a dynamical character. The spin symmetry 
of antiparticle spectrum is also well reproduced in the present calculations. 
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